active_sites = [1 2 4 7 8];
passive_sites = setdiff(1:9,active_sites);
control_sites = 10;

for site = 1:9
    for time_pt = 1:7
        cur_vel = -proj_meta(site).rd(1,time_pt).velM_smoothed;
        cur_vel = cur_vel-0.001;
        cur_vel(cur_vel < 0) = 0;
        auc_fb(site,time_pt) = sum(cur_vel(1:7500));
        auc_pb(site,time_pt) = sum(cur_vel(7501:15000));
    end
end

mean_auc_fb_act = mean(auc_fb(active_sites,:));
mean_auc_fb_pas = mean(auc_fb(passive_sites,:));
mean_auc_pb_act = mean(auc_pb(active_sites,:));
mean_auc_pb_pas = mean(auc_pb(passive_sites,:));

figure;hold on
plot(mean_auc_fb_act)
plot(mean_auc_fb_pas, 'r')
plot(mean_auc_pb_act, '--')
plot(mean_auc_pb_pas, 'r--')